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We recalculate the first analytic correction beyond PFA for a sphere in front of a plane for a 
scalar field and for the electromagnetic field. We use the method of Bordag and Nikolaev [J.Phys.A, 
41 (2008) p. 164002]. We confirm their result for Dirichlet boundary conditions whereas we find a 
different one for Robin, Neumann and conductor boundary conditions. The difference can be traced 
back to a sign error. As a result, the corrections depend on the Robin parameter. Agreement is 
found with a very recent method of derivative expansion. 
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I. INTRODUCTION 

The proximity force approximation (PFA) is one of the basic methods for calculation of Casimir and van der Waals 
forces between non-planar surfaces. Although it appeared more than 70 years ago, corrections beyond PFA, which 
gives the leading order only, were not known. Only with the new method of the Casimir force calculation, based on 
the multiple scattering formalism in conjunction with Krein's formula [lj or on a path integral quantization with a 
partial wave expansion it became possible to calculate corrections beyond PFA. Using an asymptotic expansion, 
analytic corrections were calculated, first for a cylinder in front of a plane Q, later for a sphere in front of a plane 
[1, At the same time, numerical approximations to these corrections were obtained @. 

The method, used in [J 0], results in a representation of the vacuum energy for a sphere in front of a plane (see 
Fig. 
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where £ is the imaginary frequency and M is composed of the scattering T-matrix of the sphere (see below). In (JTJ) , 
the trace is over orbital momenta, 1 = 0,1,... and m = —I, . . . , I. In case the ratio 

£ =R (2) 

of the separation to the radius of the sphere is not small, the sums and the integral in ([T]) converge rapidly allowing for 
an easy numerical evaluation of the vacuum energy and of the Casimir force. On the contrary, if e becomes small, i.e. 
for close separation, higher orders of the orbital momenta need to be accounted for. As a consequence, the numerical 
evaluation could be done down to e ~ 0.1 only. 

The analytic method makes an asymptotic expansion of ([T]) for e — > 0. This involves a substitution of the sums 
by integrals, a change of variables and an asymptotic expansion of M (involving Bcssel functions and Clcbsch-Gordan 
coefficients). 

For a cylinder in front of a plane, this method was verified independently in Q (and generalized to finite temper- 
ature). Regrettably, for a sphere in front of a plane, for Neumann boundary conditions and for the electromagnetic 
case, the results of @, Q could not be confirmed. In fact, a sign error was found and it is the aim of the present paper, 
to point to the place where it occurred and to give the correct results following from this method. 
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In the following two sections the corrections beyond PFA are recalculated for the scalar field and for the electro- 
magnetic field. The last section contains the conclusions together with a comparison with further methods. 
Throughout the paper we use units with h= c = 1. 



II. SCALAR FIELD 




FIG. 1: The configuration of a sphere in front of a plane. 

Following let the radius of the sphere be R, and the distance from the center of the sphere to the plane be L. 
Then the distance between the sphere and the plane is L — R, which we denote by d (see Fig. [T]). 
For a scalar field tp, the Casimir interaction energy is given by @, 
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EXY = ^T ^Trln(l-(-inV Y (£)). (3) 
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Here we follow the notations in Q. The first index denotes the boundary conditions on the plane. For Dirichlet 
boundary conditions, X = D and x = 0. For Neumann boundary conditions, X = N and x = 1. The second index 
denotes the boundary conditions on the sphere. Y = D, N and R for Dirichlet, Neumann and general Robin boundary 
conditions respectively. The trace Tr is the orbital momentum sum, 



oo oo 



ro=— oo l=\m\ 



and the matrix N is given by 



i+i' 



N M® = Jzk E Kl»+i/2(XL)H l u,df(SR). (4) 
' 5 i»=\i-i>\ 



Here. 



< - V(» + !)(!«' + !)(«" + 1) (o o o) (m -m o) ' ( 5 ) 

involves 3j-symbols and the function dJ(£,R) depends on the boundary conditions on the sphere. For Dirichlet 
boundary conditions, i.e., < P\ r _ R = 0, it is given by 



/2 (6R)' 
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whereas for general Robin boundary conditions with parameter a, i.e., rd r tp + aip\ r _ R = 0, it is given by 



uI l+1/a {ZR) + ZRIl +1/2 (ZR) 



uK l+1/2 m)+aRK' l+l/2 my 

The parameter u is related to the Robin parameter a by u = a — 1/2. a = or equivalently, u = —1/2, corresponds 
to the Neumann boundary conditions. In these formulas, I v {z) and K v (z) are the modified Bcsscl functions. 

The asymptotic expansion for the Casimir energy when e = d/ R <C 1 can be computed in the same way as explained 
in [H, H[. First make a substitution £ M- and expand the logarithm in © to obtain 

^ = -^- R YPk^l « £ E HE IMU- w 

s=0 1/0 m=-oo;=| m | yj = li J= | m |y \i=0 / 

The main contribution to the energy comes from I ~ l-y ~ . . . rv l s . Replacing = 1, . . . , s by I + i,, one can rewrite 
@ as 

* XY -S5 E ^rr- *EE HE n <u + ,„, - m 

s=0 1/0 Z=0m=-i V=l^ = | m |_{/ \i=0 / 

with the understanding that Iq = l s +i = 0. Using the Debye asymptotic expansions for the modified Bessel functions, 
one can show that the leading contribution comes from I ~ e _1 , U ~ i m ^ and £ ~ e _1 . This motivates to 
replace the summations by corresponding integrations, dropping only exponentially small contributions, and to make 
the substitutions, 



t i , tr tr 

m = \ — fj,, 



e e 



lAt 

-rii [i = l,. 



(8) 



One then obtains the following expression for the Casimir energy, 

where 
and 

2 Y =n(\/?^.H«-. + ,)' <«» 

For the first two leading terms of the Casimir energy, one can set e — > + directly in the integration limits for /x and 
rij, and expand Z up to terms of order e. For the term H\ v defined by J5]), it is nonzero only if I + 1' + I" is even. 
Using the substitution 

I" = 1 + 1' - 2u, 

the summation over I" in ([4]) becomes summation over v from v = to v = min{l,l'}. Since 1,1' ~ e _1 , the upper 
limit for the ^-summation can be replaced by oo. For the 3j-symbols, the small e asymptotic expansion has been 
derived in For the modified Bessel functions, the small e asymptotic expansion can be obtained using the Debye 
asymptotic expansions. In the case the sphere is imposed with Dirichlet boundary conditions, i.e., Y = D, the small 
e asymptotic expansion is given by 



-2t-{n-n'Y 



/oo 7 °° 2u / 1 \ v 
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One can then sum over v and compute the Gaussian integration over rj, which give 

-»-(„-„')»-„> U + (V2).D ^ (l),p \ 




The functions and arc given in Substituting ([IT]) into (|10[) . one can expand Z up to order s, which, 

upon substitution into © gives 



x e" 1 ' 1 ^1 + 

where r?i = X)i=o( n » ~ ™i+i) 2 , and 



i=0 




(12) 



« D = E ^WSE+E^&i- ( 13 ) 

0<i<j<s i=0 

The integrations over rii, i = l,...,s, can be performed in the same way as explained in with the help of a 
machine. The term proportional to y/e drops out since it is odd in one of the n,'s. The integrations over /x, r and t 
are straightforward. In [J], the following result was obtained for the case where Dirichlct conditions are imposed on 
the plane and the sphere (XY = DD): 



E BB = - 



16 ^E(71 1 tf( 1 + I + ---)=-t^(i + | + ---)- ( 14 ) 



Here the known formula 



°° 1 7T 4 

^ (s + l) 4 = C(4) = 90 



has been used. For the case where the plane is imposed with Neumann conditions and the sphere is imposed with 
Dirichlct conditions (XY = ND), we observe from (|12[) that the only difference is that now the summation over s is 
alternating in sign. Hence we obtain immediately 

P ND R f (-1) S+1 f, \ ^R (, \ 

s— 

where now the formula 

(-l) s+1 7 it 4 



E 



- (s + l) 4 8 90 



s=0 



is used. The results (|T4|) and (fl~5|) have been obtained in Q and Q. 

Next we consider the case where the sphere is imposed with general Robin boundary conditions. Following we 
observe that the Dcbyc asymptotic expansions give 



ul v {vz) + vzl' v (yz) =v\J\ + z 2 l v (yz) 



v zVI+^ i + f + o(^)J (lg) 

«^ ( „,( 1+ ^ + -^_, +0 (^)), 



uK v (vz) + vzK' v {yz) = - v\/l + z 2 K v {vz) I -- 



(17) 
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where 



Ml 



yrr^V8 24(1 + z 2 
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VTT^ V 8 24(1 + z 2 ) 



(18) 



We remark that in cq. (16) of which corresponds to cq. (|17[) here, there is a sign error in u which leads to results 
different from what we are going to obtain below. From (|16[) and (|17l) . we have 



1 



ul v (yz) + vzl' v {vz) 

uK v (yz) + vzK' v {yz) K v {vz) (\_ a ~"i 



K v (vz) 



Thus we obtain a dependence on the Robin parameter u which is missing in Q ■ 



To proceed, we mention that passing from Y = D to Y = R, one has a relative minus sign in front of Nf\,. The 



term dn/nLi is n °t changed, i.e., a^ 1 /^^ = ajj 1 /^^, and the term anl ni+1 is changed by 



fl/2),R _ Jl/2),D 



(1) 



a (i).R _ a W< D 

,ni-|_i ,rii+\ 



(r 2 - 1) 2w 



Since the change in the sign of Ni i> gives rise to a factor (— l) s+1 to Z, it can be compensated by the change in sign 
of the term (— V) x ( s+1 ) when one passes from X = D to X = N. Therefore we obtain from ([12")) that 



R 
4^d 2 



y^- r*te-*^ I dTV? / 



i) f 1 dry/? f°° dn ^ (s+1)/t 



n 



<7n 



x e 



E 

i=0 



(t 2 - 1) 2u 



This integral is easy to compute and we find that 



R 



y^^ 

16nd 2 ^ (s + l) 2 

s— 



e 2(6u-l) 



Therefore, 



XR 



1? 



A; V(s + i) 4 



167Td 2 ^ V(S+ l) 4 
s— x 



i 



2(6u- 1) 



7T R 

1440d 2 



1 10(6u-l) 
3 + V 2 



where we have used 



E r^T72 = c(2) = \ 

^(s + l) 2 6 



Similarly, we have 



R y> e(-l) s+1 2(6M-l) 
Wnd 2 ^ (s + l) 2 3 

s— 



This implies that 



DR. 



16nd 2 ^ \ (s + 1) 4 V 3/ 

s— 



1k a R 
11520d 2 



(-!)' 
(* + !)■ 

1 { 40(6u - 1) 
3 



e(-l) s+1 2(6u- 1) 
(s + l) 2 3 
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where the formula 

^ (s + 1) 2 ~~ 2 6 

s— 

has been applied. 

The results in this section can be summarized as 



E 



DD 



= 1 



pDD - ' o 
PFA ° 

1 + - 



PND " 1 q 
PFA 



E UR /l 80(3a-2) 



£ NR i i /'l 20(3a-2) 



^PFA \3 

Here, Epp A is the leading term which coincides with the proximity force approximation, 

ttiDD _ 77.NR _ 71 " 3 ^ pND _ pDR _ 7lT 3 R 

LpFA ~ EpFA - ^1440^' pfa ~ EpFA ~ Tl520^- 

Setting a equal to zero, which corresponds to the Neumann case, we find that 

E BN (\ 160\ 

l+{--^)e+...~l- 1.98 e- 



E™ A I 3 7^ 

£ NN (\ 40 s , 

=1+ e + ...~ 1-3.72, 



tpNN l o 2 

PFA V n 



(19) 



(20) 



We mention that the corrections involving Neumann boundary conditions on the sphere are still quite large. 



III. THE ELECTROMAGNETIC FIELD 



For electromagnetic field, if the sphere and the plane are both perfectly conducting, the Casimir energy is given by 
SI: 



1 r°° 

E EM = — d£Trln(l-N(0). 
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The trace Tr is 



Tr= E E 

m=-ao i=max{l,|m|} 

where the trace tr on the right hand side is the trace over 2 x 2-matrices - the components of N given by 



i+i' 



rJ » A\ tl , ma fdrm o 



with 



ll"(l" + 1)- 1(1 + 1) -l'(V + 1) ~ 2m£L 

A U' = o /,„ . TfTTTT, . =ff > A M' " 



^i(i + i)v(v + i) ' vf+W+i) 

and 

,TE^ _ WCg) jTM /V- d\ = ^ + l/ 2 (^)+W +1/2 (^) 

' ^ j ^ + i/2(^)' ' ^ ' ±K l+1/2 (ZR)+tRK> l+1/2 (liR) 
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Notice that dJ E and dj correspond, respectively, to df and df- with u = 1/2 in the scalar case. 

Proceeding as in the previous section, we find that the Casimir energy can be written in the same form as ©, with 
Z now given by 

z = ^fl(^^ + ia + i i+ X ( 21 ) 

To expand Z up to terms of order e, we only need, in addition to the scalar case, to expand the diagonal term Aj ( , 
up to terms of order e, and the off-diagonal term Aj up to terms of order s/e. Using the substitutions ©, we have 

j" - - ■ • • 11" 



A M , =1 + A„,„> (i/)e + . . . , \ n ,n> {v) 

A; J' =A„ n ' %/£" + ■ ■ • 7 A„ : „< : = 



frr tr' 



2/iVT 



r 



2 



Comparing to the scalar case, it is easy to see that the asymptotic expansion for Ni t i' is given by 

( /i n\ / (V 2 )>TE \ \ / (l).TE > 

e -2f-(n-n') 2 -(i 2 /T I ' 1 U > ■ l~ ' A "."' 1 , - I a -n. n' + A r: 




where 



for i = 1/2 or 1; and 



(i),TE _ (»),D (i),TM _ (i),R| 



The asterisks denote terms contributing to higher orders. 
Substituting into (|21j) and taking the trace, we find that 



\j=0 / V i=0 0<i<j'<s 



Here a TE = a D (eq. (fT3")) ) and a™ is defined in the similar way. The term of order y/e will drop out after integration 
with respect to rij. Substituting into (|9]) and compare to (fT2|) . one finds that the first two leading terms of the 
electromagnetic Casimir energy can be written as 



. first two leading 1 f first two leading 
\ermsof£: DD + terms of £ NRI 



V 2 , 



where 



^- 4 ^ 2 ^ s +iy JoVi^J-ooV^ fv—v^ 

(23) 

0<i<j<s 

Up to this point, it seems that we haven't done anything much different from the one presented in the paper [f|. 
The major difference here is that we only keep the terms in A' {/ and Aiji that will contribute terms up to order e for 
Z. In the following, we are going to see that this simplifies the computations. 

Using 



E 



VT) 2v /1-rV 1-t 2 fl-T 



i/! VI + tV 1 + r 77 6XP I 1 + r 77 

('=() 
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it is straightforward to compute \ n . n > given by ([22]) . We find that 

2(1 -r 2 ) 2 



ir 3 



1 

ir 2 ' 



Then 



i,n i+1 T f 



i=0 



0<i<j<s 



=2(.s + l) 



2(1 



ir 3 



4(s + 1) 2 (1-t 2 ) 
ir 3 ' 



2(5+1) 
tT 2 ' 



/ 2/xvT 



(24) 



■3/2 



Substituting into (|23[) and carrying out the /x-integration, it is seen that the singularity for r — > cancels between 
both contributions in the last line in (|24[) . This is equivalent to the compensation of the logarithms observed in [j| 
(see remark after eq.(135)), however now the remaining contributions come out different. We find that 



AE 



R 



X 



R 7T 



47rrf 2 ^ (s + l) 2 47rd 2 6 



■£. 



Using the results for i? DD and £' NR | u _ 1 ^ 9 from the previous section, we finally obtain the first two leading terms for 
the electromagnetic Casimir energy as 



E 



EM 



ir 3 R 
~ 1440d 2 

TT 3 R 

720d 2 
il - 1.69e 



1 + r 



TT A R 

' 1440d 2 
20' 



20 

TT-- 



R IT 

Aird 2 ~6 



■e + 



(25) 



This result agrees with that obtained recently in [8[ using derivative expansion. As has been observed in [8[ , the first 
two leading terms of the EM Casimir energy for perfect conductors turns up to be the sum of the first two leading 
terms for the DD and the NN scalar Casimir energies. 



IV. CONCLUSIONS 



The results from the calculation of the first correction beyond PFA are summarized in Eqs. (fH)]) , (f2"0"]) for the scalar 
field and in (j25j) for the electromagnetic field. These formulas substitute the final formulas Eq. (3) in Q and Eqs. 
(43) and (136) in @. It can be seen that the resulting numbers change to some extent and that the logarithmic 
contributions disappeared. 

Very recently, in [l(| and in another method for obtaining analytic corrections beyond PFA was found. It 
consists of an expansion of E, ([1]), which is perturbative in the gradients of the height profiles of the interacting 
surfaces. Intuition predicts that neglecting the derivatives, a resummation of the perturbative expansion in the height 
profiles, should reproduce PFA, while the inclusion of gradients should give the first corrections beyond PFA [10(. 
Indeed, in (Toj . the analytic result for Dirichlet boundary conditions was obtained confirming the first line in Eq. 
(fT9"|) . All other cases were computed in 0] and coincide with Eqs. (fTi?|) . (|20p and ([2~5j) . In Q, a Pade resummation 
of the asymptotic large distance expansion, matched with the first PFA correction, is performed and yields excellent 
agreement with the numerical results in 0] for DD, NN and EM boundary conditions. In addition, fits to the numerical 
data in Q, that take into account logarithmic corrections to PFA at second order (in distinction to @), yield results 
that are consistent with the results given here for the three types of boundary conditions. 
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